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Abstract 



Wc investigate conditions for the existence of the hmiting conditional distribution 
, of a bivariate random vector when one component becomes large. We revisit the 

Ph ' existing literature on the topic, and present some new sufficient conditions. We con- 

, centrate on the case where the conditioning variable belongs to the maximum domain 

of attraction of the Gumbel law, and we study geometric conditions on the joint dis- 
tribution of the vector. We show that these conditions are of a local nature and imply 
asymptotic independence when both variables belong to the domain of attraction of 
an extreme value distribution. The new model we introduce can also be useful to 
OA ' simulate bivariate random vectors with a given limiting conditional distribution. 

> 

o . 
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^' 

^ ■ 1 Introduction 

o ■ 

In many practical situations, there is a need of modeling multivariate extreme events. 
^ , Extreme means, roughly speaking, that no observations are available in the domain of 

^ I interest, and that extrapolations are needed. Multivariate extreme value theory provides 

an efficient mathematical framework to deal with these problems in the situation where 
the largest values of the variables of interest tend to occur simultaneously. This situation 
is referred to as asymptotic dependence in extreme value theory. In this case, probabilities 
of simultaneous large values of the components of the vector can be approximated and 
estimated by means of the multivariate extreme value distributions. In the opposite case of 
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asymptotic independence, the approximate probability for two or more components being 
simultaneously large given by the standard theory is zero. Refinements of the standard 
theory are thus needed. 

One refinement is the concept of hidden regular variation introduced by Resnick 2002l | . 
Another approach is studying, if it exists, the l imiting distribution of a rand om vector 
conditionally on one component being large. See iHeffernan and Resnicld 20071]. Formally 
stated in the bivariate case, this corresponds to assuming that there exist functions m, 
a and tp, and a bivariate distribution function (cdf) K on [0, oo) x (—00,00) with non 
degenerate margins such that 



lim F{X < t + '4){t)x ; Y < m{t) + a{t)y \ X > t) = K{x, y) 



(1) 



at all points of continuity of K. Das and Resnicj 20081 ] introduced the terminology of 
conditional extreme-value (CEV) model. Note that Condition ([T]) implies that X belongs 
to a max-domain of attraction. More properties can be found in Section 13. H where in 
particular the relationship between CEV models and usual multivariate extreme value 
(EV) models is explicited. Statistical applic ations of the conditional rn odel on various 
domains are discussed in several papers: see IHeffernan and Tawnl 200411 f or a study on 



air quality monitoring, lAbdous et al.l [200811 an d 



Fougeres and Soulierl " 2008 ] for an insight 



into financial contagion, and Das and Resnick 20091 ] for an application on Internet traffic 
modeling. 

An important issue that must be addressed is to study models under which Condition ([T|) 
holds. The aim of this contribution is to review existing models and exhibit new ones 
satisfying ([1]). We restrict our attention to bivariate random vectors for simplicity of 
exposition. We focus on the case where the conditioning variable belongs to the domain of 
attraction of the Gumbel distribution. The reason for that is that, as mentioned later, this 
situation is, in the models we consider, strongly related to the asymptotic independence, 
which is precisely the case where there is an advantage to work with CEV models instead 
of EV models (cf. Section 13. ip . Our work is an attempt to motivate the use of CEV 
models by exhibiting a class of bivariate models that satisfy Condition ([T]). 

In Section [21 we review the existing literature. In Section [3l we study new models for 
bivariate vectors {X,Y) with a polar representation R{u{T),v{T)) where R is a nonneg- 
ative random variable in the domain of attraction of the Gumbel law, independent of the 
random variable T, and the functions u and v are a parametrization of a certain curve. 
This model includes many models already studied in the literature, and in particular the 
bivariate elliptical distributions. Our main result (Theorem [3] in Section [3]) shows that the 
local geometric nature of this curve around the maximum of u determines the existence 
and form of the limiting distribution in ([Tj). Ou r result covers situations that are more 



general than the results of lBalkema and EmbrechtsI 2007[ ]. In particular, as a consequence 
of their local nature, our assumptions do not imply that the conditioned variable (Y) be- 
longs to the domain of attraction of a univariate extreme value distribution. Thus these 
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polar distributions may not be imbedded in a standard multivariate extreme value model. 
But when they are, we show that they are asymptotically independent. In order to prove 
this, w e extend Das and Resnick . 20081 . Proposition 4.1]. Finally, following Abdous et al 



2008l | we also study a second order correction for the asymptotic approximation ([T]). Sec- 



tion [4] contains the proof of Theorem [3l Some additional auxiliary results are given and 
proved in Section [5l 



2 Elliptical and asymptotically elliptical distributions 



Early results providing families of distribu tions that satisfy ([T]) were obtained bv lEddv and Gale 
1981 ] for spherical distributions and by Berman 19831 ] for bivariate elliptical distribu- 
tions. Mu ltivari a te elliptical dist r ibutio ns and related distributions were investigated by 



Hashorva 20061 ] . Hashorva et al. 2007 ]. One essential feature of elliptical distributions 
is that the level sets of their density are ellipses in the bivariate case, or ellipsoids in 
general. Such geom etric consider ation s have been deeply inves t igated and generalized in 
many directions by iBarbd 20031 ] and iBalkema and EmbrechtsI 20071 ] . It must be noted 
that these geometric properties will be ruined by transformation of the marginal distribu- 
tions to prescribed ones. Another specific feature of the elliptical and related models is 
that the property of asymptotic dependence or independence is related to the nature of the 
marginal distributions. If they are regularly varying, then the components are asymptoti- 
cally dependent; if the marginal distributions belong to the maximum domain of attraction 
of the Gumbel distribution, then the components are asymptotically independent. 

We start by recalling some definitions that will be used throughout the paper. A 



nondecreasing function g is said to belong to the class F or to be F- varying [Resnickl . 11987 
Definition 0.47], if there exists a positive function tp such that 



lim 



g{x + 'ipix)t) 
9{x) 



It is well known (cf. De Haan and Ferreiral 20061 . Theorem 1.2.5]) that a random variable 
X, with cdf F and upper limit of the support xi, is in the max-domain of attraction of 
the Gumbel distribution if and only if 1/(1 — F) is F-varying, i.e. 



I - Fix + iljix)t) 
lini : = e 



1 - F{x) 



(2) 



The function ip is called an auxiliary function. It is defined up to asymptotic equivalence 
and necessarily satisfies ip^x) = o{x) if xi = oo and ip{x) = o(xi — x) if xi < oo. In the 
sequel, for notational simplicity, we only consider the case xi = oo. The modifications to 
be made in the case xi < oo are straightforward, and the main change is that the rates 
of convergence must be expressed in terms of ip{x) instead of '4>{x)/x. The function ip is 
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self-neglecting (or Beurling slowly varying, cf. Bingham et al. 1989I . Section 2.11]), i.e. 
for aU t > 0, 



lim 



ip{x) 



1 . 



Two random variables X and Y such that {X, Y) belongs to the maximum domain 
of attraction of a bivariate extreme value distribution G are said to be asymptot ically 
independent if G has independent marginals. (See e.g. jPe Haan and Ferreiral . |2006| . Sec- 
tion 6.2]). 

In the following two subsections, we recall the results on elliptic dis t ribut ions and we 
state a bivariate version of a general result of Balkema and Embrecht j 2007 ] which pro- 
vides geometric sufficient conditions for ([T]) to hold. We also point out and illustrate the 
local nature of the sufficient condition formulated by iBalkema and EmbrechtsI 20071 ] , in 
the case of F-varying upper tails. 



2.1 Elliptical distributions 

Consider a bivariate elliptical random vector, i.e. a random vector {X, Y) that can be 
expressed as 



iX,Y) = ii(cose,/9cose + crsine) 



(3) 



with = 1 — p^, in terms of a positive random variable R called "radial component" and 
an "angular" random variable G uniformly distributed on [0, 27r). The following result 
wa ,s originally proved in the case p = as a technical lemma unde r restrictive con ditions 
Eddy and Galel 1981]. The gene ral r esult was first prove d in Berman 1983] in the 



m 



BermanI |l992l ] and lAbdous et al.l \200^ ) and 



Hashorva 



20061 ] 



m a 



bivariate case (see also 
multivariate setting. Throughout the paper $ will denote the cdf of the standard normal 
distribution. 



Theorem 1 dBermanl \l98^ ). Let iX,Y) be an elliptical random vector as defined in 
1^. // the radial component R is in the domain of attraction of the Gumbel law, i.e. its 
survival function 1 — H satisfies ^ with auxiliary function tjj, then X and Y are also in 
the domain of attraction of the Gumbel law and 



lim P(X < t + ilj{t)x ,Y <pt + ^/t'4){t)y \ X > t) = {I 



'My) 



(4) 



Comments on Theorem 1 The fact that X has F-varying upper tails follows from 
by taking y = oo. Since the same result holds with reversed roles for X and Y, the 



consistency result jPas and Resnickl . l2008l . Theorem 2.2] implies that (X, Y) belongs to 
the domain of attraction of a bivariate extreme value distribution. Moreover, since the 



4 



limiting distribution in (jl]) has two independent marginals, X and Y are asymptotically 
independent, i.e. the limiting extrern e value distribution is the product of its marginal 
distributions. See also Hashorva 20051 . Section 3.2] for a proof of this property in a related 
context. In the case where R has a regularly varying tail, the limiting distribution is not a 
product and the vector {X,Y) is asymptotically dependent. See e.g. [Abdous et al.l . l2005l . 
Theorem 1, part (i)]. As mentioned in the introduction, we do not develop this case. 

If the radial component R has a density /i, then the vector [X, Y) has the density / 
defined by 



f{x,y) 



h{^x^ + {y - px^/a'^) 
y^x^ + (y — pxY ja'^ 



The level lines of the density are homothetic ellipses + (y — px)"^ /a"^ = (? 




Figure 1: Level lines of the density of an elliptical distribution. The slope of the straight 
line is /9 = .6. 

This result can be generalized in two directions: either by weakening the assumptions on 
the level lines of the density or by extending the representation ([3|) . The first generalization 
will be considered in the following subsection, and the second one in Section [3l 



2.2 Asymptotically elliptical distributions 



In this section, we state a bivariate version of iBalkema and Embrechtd 120071. Theorein 
11.2]. We first need the following definition taken from iBalkema and EmbrechtsI 20071 . 
Section 11.2]. 
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Definition 1. A function L :M? ^ R+ belongs to the class C if for all {x,y) E M^, 
for any norm \ \ • \ \ on . 

Assumption 1. The random vector {X,Y) has a density f such that 

f{x,y)=e-'^^'y^L{x,y), (5) 

where L £ C, and the function I satisfies: 

I{x,y) =pon{x,y) , 
^ ^ r ds 
^^^^ = JoW)^ 

ip is absolutely continuous with lima;^oo = and n : ^ M is 1- homogeneous, r? 
is twice differentiable and the Hessian matrix of n is positive definite. 

Theorem 2 ( Balkema and Embrechti 2007]). Under AssumptionUl X and Y are in the 
domain of attraction of the Gumbel law and satisfy ^ with auxiliary function ijj, are 
asymptotically independent, and there exist real numbers p and a such that 



lim P(X < t + Tp{t)x ]Y <pt + ay/tiHt)y \ X > t) = {I - e-^)$(y) . (7) 



This result shows that Assumption [T] is a sufficient condition for the hmit ([T]) to hold, 
with m{x) = px, a{x) = a^xil){x) and K{x,y) = (1 — e~^')<I>(y). The constants p and a 
are characterized by the second order expansion of the function n 

n{l + x,p + ay) = l+x + y'^/2 + o{x^ + y'^) . 

This condition implies that the tangent at the point (l,p) to the curve ?^(C,C) = 1 is 
vertical. Note that the level lines of the function n are not those of the density / defined 
in ([5]), unless the function L is constant, but, loosely speaking, the level lines of / converge 
to those of n. Theorem [2] implies Theorem [T] when the radial distribution is absolutely 
continuous. 

Example 1. Let h and g be density functions defined on [0, oo) and [— 7r/2, 7r/2], respec- 
tively. The function / defined by 

fi^.y) = ^^/V^^'^'^^'Ju'^y, 9 ° arctan((y - px)/x./Y^^) (8) 
+ {y- pxY/{l - p^) 

is then a bivariate density function on M^. If is a constant, then / is the density of an 
elliptical vector. If h can be expressed as in ([6]) and if g is continuous and bounded above 
and away from zero, then / satisfies Assumption [TJ Figure [2] shows the level lines of such 
a density, with p = .6, h{t) = exp(tV2)/\/^ and g{t) = c{l + [t^ - {-K/Aff}. The level 
lines seem to be asymptotically homothetic. 
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Figure 2: Level lines of the density given by ([8]). 

It is important to note that under Assumption [1] the normalizing functions m and a 
satisfy a(x) = o(m(x)), since in the present context m{x) = px and a{x) = ay/ayip{x) with 
ip{x) = o{x). This implies that only the local behaviour of the curve n{^,Q = 1 around 
the point (1, p) matters. In other words, the limit ([7]) still holds if (X, Y) is conditioned to 
remain in the cone {{p — e)x <y< {p + e)x} for any arbitrarily small e > 0. This suggests 
that Assumption [1] must only be checked locally to obtain the limit ([7|). 

Example 2 (Mixture of two bivariate Gaussian vectors). Let B he a. Bernoulli random 
variable such that P(i? = 1) = p G (0, 1). Let X and Z be two i.i.d standard gaussian 
random variables, p ^ t £ [—1,1] and define Y by 

Y = B{pX + v^l - P^Z) + (1 - B){tX + r^Z) . (9) 

Then y is a standard Gaussian variable, and {X, Y) is a mixture of two Gaussian vectors. 
Figure [3] shows the level curves of the density function of the pair {X, Y) with p = .4, 
p = .8 and r = —.4. The density function of {X,Y) does not satisfy Assumption [H and 
Theorem [2] cannot be applied. Indeed, applying Theorem [T] to each component of the 
mixture yields 

^{Y - px< \l\ - p^z \ X>x) 

= pF{p{X -x) + Vl - p'^Z < Vl - p'^z \ X >x) 

+ (1 - p)¥{t{X - x) + (t - p)x + \/l - T^Z < y/l - p'^z \ X>x) 

~p$(z) + (l-p)l|,<p} . 
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Thus the hmiting distribution is degenerate, with a positive mass either at — oo or +00. 
However, a proper hmiting distribution can be obtained for (X, Y) conditioned to remain 




Figure 3: Level hues of the density function of the pair (X,Y) defined in Example [2j 

in C = {{x,y) G | cix < y < C2x} for ci < C2 such that p G [ci,C2] and r ^ [ci,C2] 
Denote P(- | {X,Y) G C) by Pc(-)- Then 



lim ¥c(Y < px+ a/1 - p^z \ X > x) = . 

To prove this claim, we assume without loss of generality that p = 0. Then 

F{Y <z; {X,Y) gC \ X > x) 



¥c{Y <z\X >x) 
For fixed z and x > z/c2, it holds that 



p((x,y) £C\x >x) 



F{Y <z; {X,Y)£C\X >x) 
= p¥{ciX < Z < z\X >x) 



+ (1 - p)F{ciX <tX + \/l -t2Z < z I X > x) ~ p^{z) . 
Since p G [ci, C2] and r ^ [ci, C2], and since 2; < C2X, it is easily obtained that 

lim ¥{ciX <Z<z\X>x) = ^{z) , 



lim P(ciX < rX + \/1-t'^Z <z\X> x)={) . 
Thus lim^^ooIP(i^ <z; {X,Y) eC \ X > x) =p^iz), which proves 1^. 



(10) 



3 Bivariate vectors with polar representation 



In this section, we show that Theorem [T] can be extended from eUiptical distributions 
to more general bivariate distributions that admit a radial representation R{u{T),v{T)) 
where R and T are independent, T is not uniformly distributed and u and v are more 
general functions than in the elliptical case. We start by collecting the assumptions that 
will be needed. 

Assumption 2. 

A The function u : [0, 1] [0, 1] is continuous, has a unique maximum 1 at a point 
to G (0, 1) and has an expansion 

u{to + t) = I - £{t) (11) 

where i is decreasing from [— e, 0] to [0, r/_] and increasing from [0, e] to [0, rj^] for 
some e,ry_,ry_|. > 0, and regularly varying at zero with index k > 0. The functions 
■ [0,1]'-] ~^ and l"^ : [0,r7+] [0, e] respectively defined as i^{s) = sup{t E 

[— e, 0] : l{t) < s} and ^+ (s) = inf{t E [0, e] : l{t) > s} are absolutely continuous and 
their derivatives (£!r)'(s) and are regularly varying at zero with index 1/k— 1. 

B The function v defined on [0, 1] is strictly increasing in a neighborhood of to, u(to) = P; 
and the function t i— > v{to + t) — p is regularly varying with index (5 > 0. Its inverse 
is absolutely continuous and its derivative is regularly varying at zero with index 
1/5 -I. 

Assumption 3. The density function g : [0,1] is regularly varying at to with index 

T > — 1 and bounded on the compact subsets of [0, 1] \ {to}- 

Theorem 3. Let R be in the domain of attraction of the Gumbel law with auxiliary 
function ip, i.e. its distribution function H satisfies l^j- Let T be a random variable 
that admits a density g that satisfies Assumption Let the functions u and v satisfy 
AssumptionlMwith 6 < K. Define {X,Y) = R(u{T),v{T)). Then, 

(i) the random variable X is in the domain of attraction of the Gumbel law and there 
exists a function k regularly varying at zero with index (1 + t)/k such that 

F{X > x) k{tl;{x)/x)H{x) ; (12) 
(a) there exists a function h regularly varying at zero with index 6/k such that for all 

yem, 

lim F{X < t + ^p{t)x ,Y < pt + th{i^{t) /t)y \ X > t) 

= (l-e-^)i/./5,(i+.)/5(y), (13) 
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with 

J — OO I I 

The proof of this result is in Section HI One of its main ingredient is the fact that the 
tail of R is F-varying. This implies that the normalizing function a{x) = xh{ip[x) / x) is 
o{x). As a consequence, only the local behavior of u and v around (l,p) matters. This is 
similar to what was observed under Assumption [TJ 



Comments on Theorem [3] 



(i) The case u{t) = cos(t) and v{t) = sin(t) was considered bv iHashorval |2008l |. The 
difference between the present results and this reference is not only that we consider 
more general functions u and v, but more importantly that we point out the local 
nature of the assumptions on u and v. 

(ii) Theorem [3] handles situations where the assumptions of Theorem [2] do not hold. In 
some cases, the limiting distribution is nevertheless the Gaussian distribution and the 
normalization is the same as in Theorem[3l see Example[3l In other cases, the limiting 
distribution and the normalization differ from those that appear in Theorem[2l There 
are two reasons for this: the density of T can vanish or be unbounded at zero, or the 
curvature of the line parameterized by the functions u and v at the point (li(to), ^^(to)) 
can be infinite or zero. This is illustrated in Example HI 

The asymptotic distribution -^^/^^(i+r)/^ is of the so-called Weibull type with shape 
parameter k/5. This ratio characterizes the geometric nature of the curve t i— > 
{u{t),v{t)) around t^ and is independent of the particular choice of the parametriza- 
tion t I— > [u{t)^v{t)). Its right tail is lighter than the exponential distribution. The 
behavior of the density of T has influence only on the less important parameter Q. 
The normalizing function h also depends only on u and v. 

(iii) Denote a{x) the normalizing function in (jl3j) : a{x) = xh{ip{x)/x). The assumption 
5 < K implies that il){x) = o(a(x)), so that (AT — x)/a{x) converges weakly to zero 
given that A" > x as a; — > oo. Hence 



Y — pX 

lim IP ( ^^^^ < y I A > X ) = H^/s^^i+r)/s{y) 



This implies in particular that the case p ^ can be deduced from the case p = 
by a linear transformation. 

Since the function h is regularly varying with index 6/k > and since ip is self- 
neglecting, it also holds that a{X)/a{x) converges weakly to 1 given that A > x as 
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oo, thus the conditional convergence also holds with random normalization: 



This is a particular case of iHeffernan and Resnicld [2007], Proposition 5] 



(iv) Note that Theorem [3] states that X is in the domain of attraction of the Gumbel 
law, with the same auxiliary function ijj as R. However, little can be said about the 
tails of Y other than P(y > y) < IP(-R > y/v*), where v* = max^gjg,!] ^(0- This 
is because nothing is assumed about of the behaviour of v and g around the point 
where v has a maximum. This would be needed to obtain an asymptotic form for 
the tail of Y similar to (112p . It is in contrast with the situation of Theorems [T] and O 
Nevertheless, it can first be proved that 

bvit) ~ v*b{t) (14) 

where b and by are the inverse functions of 1/P(i? > •) and 1/¥{Y > •), respectively 
(see a proof in Section [5]). This implies that if Y does belong to the domain of 
attraction of an extreme value distribution, then this distribution is necessarily the 
Gumbel law (since Y has a lighter tail than v*R and unbounded support) and {X, Y) 
belongs to the domain of attraction of a bivariate extreme value distribution with 
independent marginals, i.e. X and Y are asymptotically independent. This is shown 
in Corollary [5] below. 



3.1 Relations with CEV and EV models 



As mentioned in the introduction. Das and Resnick 2008l | referred as CEV models the fam- 
ilies o f distributions satisfying Condition ([T|). An important finding of lHeffernan and Resnick 



2007l | is that in such a model, it is not possible to transform X and Y to prescribed 



marginals Fi{X) and ^2(^)5 for given univariate cdf Fi,F2, when the limiting distribu- 
tion K is the product of its marginals. Theorem [3] provides random vectors {X, Y) for 
which the limiting distribution in ([1]) is precisely a product, so that nonlinear transforma- 
tions of the marginals to prescribed marginals are impossible. This is in contrast with the 
usual practice of standard multivariate extreme value theory. 



As noted bv iDas and Resnickl 2003 . Section 1.2], another advantage of the CEV ap- 
proach is that it does not require the assumption that all components of the vector belong 
to the domain of attraction of a univariate extreme value distri bution. The relat i onship 
between CEV models and EV models has been investigated in Das and Resnick 20081 ]. 
They proved in particular that if (X, Y) belongs to the domain of attraction of a bivari- 
ate extreme value distribution with asymptotic dependence, the n the CEV model does 
not provide any more information than the EV model. See also iDas and Resnickl 20091 . 
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Section 1.2]. Conversely. iDas and Resnickl [20091 . Proposition 4.1] gives conditions under 
which, if (X, Y) satisfies ([1]) and Y belongs to the domain of attraction of a univariate 
extreme value distribution, then {X, Y) belongs to the domain of attraction of a bivariate 
extreme value. 

The next result elucidates the relationshi p between the condition al limit of Theorem [3] 
and extreme value theory. It is similar to Das and Resnick 20081 . Proposition 4.1] but 
covers cases ruled out by this reference, as shown afterwards in Corollary [5j In the 
following, the extreme value distribution with index 7 is denoted by G^, that is to 
say G.y{y) = exp{ — (1 + jy)~^^'^} for each y such that 1 + 7y > for 7 7^ 0, and 
Go{y) = exp{-e-^}. 

Proposition 4. Assume that the vector {X, Y) satisfies HJ) and that Y belongs to the 
domain of attraction of an extreme value distribution G^ with auxiliary function ipy, i-S- 



lim P(y >t + tpY{t)y \ Y >t) 

t— »oo 



-logG^(y) . 

Assume moreover that for all y G M, 

6y(t) - m o + V'y o bY{t)y 



lim 



a o hx{t) 



+00 



(15) 



where hx and by are the inverse functions o/l/P(X > •), of 1/P(y > •), respectively. 
Then, (X, Y) belongs to the domain of attraction of a bivariate extreme value distribution 
with independent marginals, i.e. X and Y are asymptotically independent. 



Proof. To prove asymptotic independence, we must show that 

lim t¥{X >bx{t) + tlJo bx{t)x , Y > by {t) + i/jy ° by{t)y) 







(16) 



where -i/; is as in ([T]) . For any ?/ G M, define y by bY{t)+'ipYoby{t)y = mobxit) + aobx{t)y, 
i.e. 



y 



6y(t) — m o bxit) + tpy o by{t)y 



a o bx{t) 

Thus, by (fT5]) . for any z G M, we have, 

lim sup tP(X >hx{t)+^o bx{t)x , Y > by{t) +^y o by{t)y) 



t—too 



lim sup tF{X > bxit) +ipo bx{t)x , Y > mobx{t) + ao bx{t)y) 



t^oo 



< limsuptP(X > bxit) + ^ o bx{t)x , Y > m o bx{t) + a o bx{t)z) 



= ¥{X* > X , Y* > z) , 

where {X* ,Y*) is a random vector with distribution K. This can be made arbitrarily 
small by choosing z large enough, so (fT6]) holds. □ 
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A simple example is provided by the elliptical distributions with identical margins, for 
which one has: bx = by, iprix) = o{x),tpY{x) = o{a{x)) and m(x) = px, so that (fT5]) 
holds. A more general result is provided by the following corollary. 

Corollary 5. Under the assumptions of Theorem\^ if moreover Y belongs to a domain of 
attraction, then it is the Gumbel law, and il5\) holds, so that X and Y are asymptotically 
independent. 



Proof. Under the assumptions of Theorem [3l the normalizing functions in ([T]) satisfy 
m(t) = pt, il){t) = o{a{t)) and a{t) = o(t). Since R has unbounded support, then so 
has Y and since Y has lighter tails than v*R, the max-domain of attraction of Y can 
only be the Gumbel law. Thus it also holds that ipvit) = o{t). By (jl4p we know that 
6y(i) ~ v*b{t) and (fT2]) implies that bx{t) ~ b{t). Thus we have 

byit) -mobx{t) + ipY ° bY{t)y {v* - p)b{t) 



a o bxit) 

The assumptions on v imply that v* > p, so 

t-+oo aobx{t) 

Thus (HI]) holds for all y G R. 



a o bx{t) 



oo . 



□ 



Remark iDas and Resnicld 2008 . Proposition 4.1] show the same result under a condition 
which can be expressed with the present notation as lim^^oo i^Y o 6y (t)/a o bx{t) € (0, oo]. 
Under the conditions of Theorem [31 if moreover v and g satisfy some smoothness as- 
sumptions around the maximum of v similar to Assumptions l2IAIand[3l it can be shown 
that tpY ~ ip and thus lim^^oo i^Y ° byit) /a o bx{t) = 0, so that iDas and Resnicld 20081 . 
Proposition 4.1] cannot be applied here. 



3.2 Some applications of Theorem [3] 

We now give some examples of applications of Theorem [3l 

Example 3. If the density g of the variable T has a positive limit at tQ, then r = 0. If n 
is twice differentiable with u'{to) = and n"(to) / 0, then k. = 2 and if v'{0) > 0, then 
6 = 1. If these three conditions hold, the limiting distribution is the standard Gaussian 
and the normalization is a^Jxil){x). Figure [2] also illustrates this case. As will be shown 
in Section 13.31 this is actually a particular case of Theorem [2j 



Example 4. iHashorva et al.l [20071 ] have introduced a generalisation of the elliptical distri- 



butions, which they called Lp-Dirichlet distributions for all p > 0. We consider here the 
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case p > 1. Instead of being ellipses, the level lines of the density of these distributions 
have the following equation: 

with p G (—1, !)• See Figure [H To simplify the discussion, we consider the case p = 0. 
An admissible parametrization is given by u{t) = (1 - tP)^^P and v{t) = t, which yields 
6 = 1 and k = p. Thus, Assumption [1] does not hold except if p = 2, which is the elliptical 
context. If the density of T has a positive limit at zero, then r = and the cdf of the 
limiting distribution is then 

J —oo 



2pi/p-ir(i/p) 



Figure 4: The curve given by Equation (jl7p with p = 3/2, p = (left) and p = .7 (right). 

As a last remark in this section, note that the result (jl2p on the upper tail of X is 
actually a particular case of a more general result on the tail of the product of a random 
variable X in the domain of attraction of the Gumbel law and a bounded random variable 
U. Similar results exist for heavy- tailed or s ubexp one ntial distributions, for instance th e 



celebrated Breiman's Lemma (see 



Breiman 196 



xpone ntiai aistriputions, tor instance ti n 
^ and Cline and Samorodnitsky 1994l |) 



In the present context, Hashorva 20081 . Theorem 3] states a version of this result with 
U = cosG and has a density on (— vr, vr). We state it under slightly more general 
assumptions which highlight the fact that only the behaviour of U near its maximum 
must be specified. 

Proposition 6. Let R be a nonnegative random variable whose cdf H is in the domain 
of attraction of the Gumbel law. Let U be a nonnegative random variable, independent of 
R, such that U < b < oo a.s. and that admits a density g in a neighborhood of b which is 
regularly varying at b with index r > — 1. Then RU satisfies (0j with auxiliary function 
ij){x/b) and 



F{RU > x) ~ b'^T{T + 1) 



i;ix/b) 



g{{l/b + ip{x/b)/x}-^)H{x/b) 



(18) 
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Proof. The proof is along the hnes of the proof of (I12p , and makes use of Lemma El □ 



3.3 Relation between Theorems [2] and [5] 

Let X, Y be random variables whose joint density / can be expressed as 

f{x,y)=gon{x,y), (19) 

where g' is a nonnegative function on M+ such that rg{r) dr < oo, n is a positively 
homogeneous function with index 1 and the level line n{x, y) = 1 admits the parametriza- 
tion t — > {u{t),v{t)), t G [0, 1]. The change of variable x = ru{t), y = rv{t) yields, for any 
bounded measurable function ip: 

roc rl 

E[ip{X,Y)]= / ip{ru{t),rv{t))rg{r)\u{t)v'{t) - u'{t)v{t)\drdt . 
Jo Jo 

Hence {X, Y) = R{u(T), v{T)), where R has density c{g)~^rg{r), with c{g) = rg{r) dr 
and T has density c(n, v)~^\uv' — u'v\ with c(n, v) = /J \u'{t)v{t) — u{t)v' (t)\dt. 

Conversely, if {X, Y) = R{u{T), v{T)) where R and T are independent, R has a density 
h on [0, oo), and if there exists a point Iq such that u'^Iq) = 0, u"(to) > and t^'(to) > 0, 
then the function v/u is invertible on an interval around to- Let (p be its inverse, and 
define n{x,y) = x/{uo (p(y/x)}. Without loss of generality, we can assume that it(to) = 1 
and denote w(to) = P- Then n is a positively homogeneous function with index 1, and if 
T admits a density g, then {X, Y) has a density / in a cone C = {{x,y) \ (p — e)x < y < 
(p + e)x} around the line y = px, defined by 

h{n{x,y)) 
n{x,y) 

with 

g o (t){y/x) 



L{x,y) 



u'v — uv'\ o (j)[y/x) 

If the density g and the Jacobian \u'v — uv'\ are both positive and continuous at to, then 
the function L belongs to the class C. See Lemma [9] for a proof. Thus Assumption [1] 
locally holds, and Theorem [2] implies Theorem [3] in this context. 



3.4 Second order correction 



As illustrated in lAbdous et al. 2008l |. it is useful for statistical purposes to have a second 
order correction to the asymptotic approximation (|13p provided by Theorem [3l In order 
to obtain such a refinement, rates of convergence in all the approximations used to prove 
TheoremOare needed. To simplify the discussion, we will consider the following additional 
assumptions. 
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The random variable T is uniformly distributed over [0, 1]. 
There exist A > and o" > such that 



uito + t) = l 



+ o(t3), v{to + t) = p + \t + 0{t^) . 



There exist functions x ^.nd B such that 

H{x + t%lj{x)) 



H{x) 



< x{x)B{t) 



(20) 



for all t > and x large enough, where lima;^oo x{^) = 0) ^-^id B is bounded on the 
compact subsets of [0,oo) and integrable over [0,oo). 



The bound (|20p is a nonuniform rate of convergence. See lAbdous et al.l 20081 . section 2.2] 
for examples. Under these assumptions, it is possible to obtain a rat e of convergence and 
a second order correction in Theorem [3l Proceeding as in the proof of Abdous et al. 20081 . 
Theorem 3] yields 



P(F < px + -Jx%l;{x)z \ X> x) 
a 



+ 0(x(x)) + o (TVK^ 



Replacing zhy z + p\ ^ y^Tp{x)/x yields 
^(y < px+ ^y/xip{x)z + ^ipix) \ X >x^ 



^(z) + 0{x{x)) + o (VV'(x)/x) . 



This second order correction is meaningful only if xi^) = o(y^ •iIj{x)). The improvement 
is only moderate here: the bound is o{^J'4j{x)/x) instead of 0{^J'4j{x)/x)^ because we 
assumed only that n(to + t) = 1 — a'^t^ /2 + o{t^). If the expansion of u around to is 
u{tQ + t) = 1 — cj^t^/2 + O(t^), then the bound becomes 0{tlj{x) /x). This is the case for 
bivariate elliptical distributions. 
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4 Proof of Theorem [3] 



Let e be defined as in Assumption [2j Since u has its maximum at to, there exists r/ > 
such that for all t ^ [Lq — e, Iq + e], it holds that u{t) <! — ?/. Then, 

P(X>.;F>!,)=y||'ij(^vJ^)»(t)dt 

+ / H (^\/ ^] g{t)dt . 

Let r{x) denote the last term. The bound (j22p in Lemma [7] yields that for any p > 0, 

r{x) < H{x/{1 -r])) = {{tl;{x)/x}PH{x)) . 

This will prove that r{x) is negligible with respect to the first integral for which we now 
give an asymptotic equivalent. By assumption, we can choose e such that the function 
v/u is continuous and increasing on [to — e,to + e], because (5 < k, so that v/u ^ v in a 
neighborhood of to- 

If y can be expressed as y = px + o(x), then for large x, it holds that y/x £ [{v/u){tQ — 
e),{v /u){tQ + e)]. If y/x > p = {v/u){to), then there exists ti £ [to, to + e] such that 
{v/u){ti) = y/x. Thus, 

Let / and J denote the last two integrals, respectively. The successive changes of variables 
s = l/u{t) and s = 1 + uj^{x)/x yield 

J= / Hixs) ^ ' ^ g{u^{l/s))ds 

h/u{ti) s 

rx{l/u{to+e)-l}/4){x) 



ip{x) 



x 



H{x + ip{x)uj) 

x{l/u{t-,)~l}/i^{x) 

{u ^)'{l/{l+uj^{x )/x}) 

{1 + LOtjj{x) / xY' 



X ^ (\ \ .^i.r^\/Jx2 9[u (1/{1 + uJij{x)/x})) dw 



Let ki denote the function defined by 

hiio) = -c^(n^)'(l/{l + Lo})g{u- + u;})) . 
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Assumptions [2] and [3] imply that ki is regularly varying at zero with index (1 + t)/k and 
Lemma [8] yields 

/•oo 

J ki{i){x)/x}H{x) / e"¥^+^)/'^~^ dt , 

J 20 

where zq = x{l/u{ti) — l}/ip[x) and if y is chosen in such a way that zq has a finite limit 
when X — > oo. Set = /O + Then, by definition of ti, zq = x{l/u o {v/u)^{p + 
^) — l}/^{x). Assumption [2] implies that the function uo {v/u)^{p + ^) is regularly 
varying at zero with index k/5. Define an increasing function h which is regularly varying 
at zero with index (5/k by 

h{x) = (v/u) o n^(l/{l + x}) — p = (1 + x)v o + x}) — p . 

For z>0, set y/x = p + h{'il){x) / x)z. Then 

= -^h^{h{^{x)/x)z) ~ z^/' , 

POO 

J ~ ki{iP{x)/x}H{x) / e-*t(^+^)/^"^ dt . 

We next deal with the integral /, still in the case y/x > p. Noting that {y/v{ti) — 
x}/7p{x) = Zq, the changes of variables s = l/v{t) and s = {x + '4j{x)uj)} /y yield 



1= / ''^ g{v^{l/s))ds 

Jl/vih) S 

yil^{x) 



{y/v(to-t)-x}/ii)(x) 

H{x + '4j{x)uj) 
{v^)'{{y/x)/{l+ujij{x)/x]) 



{l + uji^{x)/xY 
The choice y/x = p + h{ip{x)/x)z also yields 

y/x p + h{il){x) / x)z 



9{v^ {{y / x) / {I +u:^{x)/x]))duj . 



p + h{7p{x)/x)z . 



l + uJ^l^{x)/x 1 + (jOip{x)/x 
Let the function k2 be defined by 

k2{o = ({vn'ip + m)9 ov-{p+ HO) . 

Assumptions [2] and [3] imply that /c2 is regularly varying at zero with index (l + r)/K + l — 
6/k>{1 + t)/k. Lemma El yields 

/•oo 

/ - pz^^+^'^^^'^k2{'il^{x)/x)H{x) / e-^dt = o{J) . 

The case y/x < p can be dealt with similarly and is omitted. The remaining of the proof 
of assertions (i) and (ii) is straightforward. □ 
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5 Lemmas 



The followi ng Lemma is a straigh tforward consequence of the representation theorem for 
th e class F Bingham et al. . 1989I . Theorem 3.10.8]. The argument was used in the proof 
of lAbdous et al.l 20051 . Theorem 1]. We briefly recall the main lines of the proof for the 
sake of completeness. 

Lemma 7. Let H be a cdf in the domain of attraction of the Gumhel law infinite right 
endpoint. For any p > 0, there exists a constant C such that for all x large enough, and 
all t > 0, 

H(x -\- 'ib(x)t) ^, ^ „ . ^ 

V ^H^y ) ) ^ c{l+t)-P . (21) 



H{x) 
H{ax) 



< C{iP{x)/xf . (22) 



H{x) 

Proof. If 7 = 0, the function H can be expressed as 

where lima;^oo c(3;) = c € (0, oo) and lima;^oo '0'(^) = 0- Thus, for any e > and x large 
enough, there exists a constant C such that 

c{x + il){x)t) ^ ^ tpjx + 'ip{x)t) ^^j^^^ 

c{x) " ' i^ix) ~ 

Hence 

H{x + il){x)t) 



H{x) 



This proves (j2ip . The bound ()22p follows trivially from (|2ip by choosing e <l/p and by 
setting t = {a — 1)x/iIj{x). □ 

Lemma 8. Let H be a cdf in the domain of attraction of the Gumbel law with infinite right 
endpoint. Let g be a function regularly varying at zero with index r > — 1 and bounded on 
compact subsets of (0, oo] . Then 

' H{x) 9{4,(x)lx) A 



locally uniformly with respect to z > 0. 
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Proof. Denote x{^) = ipix)/x; then lim^^oo x{x) = 0. By assumption, H{x+ip{x)t)/ H{x) 
converges to e~* and g{x{x)t) / g{x{x)) converges to , and both convergences are uniform 
on compact sets of (0, oo). It is thus sufficient to prove that 

P H{x + i,{x)t)g{x{x)t) 

hm hm sup / . , dt = , (23) 

A^co ^.^oo Ja H{x) g[x[x)) 

, , n H{x + Mx)t) gix{x)t) , , 

Urn hm sup / ^ ' ' di = o . 24 

Let e > be such that 1/e — 1 > r. By Lemma[71 for large enough x, 



/■oo 



Ax{a;) 

Since (7 is locally bounded 011 (0,oo] and —1/e + r < —1, Karamata's Theorem (cf. for 
instance 
such that 



Bingham et al. 19891 . Proposition 1.5.10]) implies that there exists a constant C 

t-^''g{t)dt < C'ixix)A)'~'/^g{x{x)A) 



r 

JA: 



lAx(x) 

P H{x + V^(x)t) g{x{x)t) ^ g{x{x)A) 

limsup / -—-^Q.t<CCA ' limsup — , , 

x^oo Ja H{x) g[x[x)) x^oo g[x[x)) 

= CC'A^-^I^^^ 

as A tends to infinity, because 1— 1/e+r < 0. This proves (p3]) . Since H{x^'4){x)i) j H {x) < 
1 and by Karamata's Theorem, we get, for some constant C, 

^ H(x + ^(x)t) , , , , , r . . ^ ^ , 

^—^-^^^g(x{x)t)dt< / g{x{x)t)dt 
H{x) Jo 

rxix)v 

= X(.x)~^ / g{t) dt < Crig{x{x)ri) . 
Jo 

Hence 

li^ li„ /' g(^f^Wt).(xWt)^, ^ _ „ 

77^0 ^^00 Jo H{x) g{x{x)) v^o 



because 1 + r > 0, which proves □ 

Lemma 9. Let i be a continuous function defined on [0,oo), bounded above and away 
from zero and with a finite limit at infinity. Define L on x by L{x,y) = t{y/x). 
Then L belongs to the class C. 
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Proof. Since £ is bounded, it suffices to prove that if the limit lim^^^^Q_,^£{{x + 0/(2/ + 
C))/£{C/C) exists, then it is equal to 1. Since moreover £ is continuous and bounded away 
from zero, it is enough to consider subsequences and to show that if ||(Cn)Cn)ll ~^ oo and 
if lim^^oo Cn and hm„_»oo Cn both exist, then lim„^oo(a; + (,n)/iy + Cn) = lim 
Three cases arise. 



(i) If lim„_oo Cn = lim„^oo Cn = oo, then 

X + Cn Cn l + x/Cn 



Cn/ Cn 



y + Cn Cn 1 + y/Cn 

(ii) If lim„^oo Cn < oo and hm^^oo Cn = oo, then 

lim (x + Cn)/ (y + Cn) = = lim Cn/Cn ■ 
n— »oo n^oo 

(iii) If lim„_^oo = oo and hm„^oo Cn < oo, then 

lim {x + Cn)/{y + Cn) = oo = lim Cn/Cn • 



□ 



Proof of igp. By assumption, F(Y > y) < H{y/v*), thus byit) < v*b{t). Denote V = 
v{T). Two cases only are possible: (i) v* is an isolated point of the support of the 
distribution of V and f{V = v*) > 0; (ii) P(y = v*) = and for any e > 0, there exists 
V e {v* - e, V*) such that ¥{V > v) > 0. 

(i) If ¥{V = V*) > and v* is an isolated point, then there exists v** < v* such that 

P(y > y) = H{y/v*)F{V = v*) + F{RV >y , V < v**) . 

Note that F{RV >y , V < v**) < H{y/v**), and since H is T- varying, H{y/v**) = 
o{H{y/v*)), thus P(y > y) ~ eH{y/v*) with e = ¥{V = v*). Since h is slowly 
varying, this implies that fey(t) ~ v*h{et) and finally that 6y(t) ~ v*b[t) as t — > oo. 

(ii) In the second case, fix some v < v* and denote e = F{y > v) > by assumption. 
Then 

P(y >y)= F{RV >y , V > v) + F{RV >y , V < v) > eF{R > y/v). 
Thus, 6y(i) > vh{et) and since h is slowly varying, 

hm mf — —— > V . 

t^QO b{t) 

Since v can be chosen arbitrarily close to v* and since we already know that 6y(i) < 
v*b{t), we conclude that limt-_»oo ^y(i)/^(^) — ■ 

□ 
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